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Abstract. Let mi, . . . , m s be positive integers. Consider the sequence 
defined by multinomial coefficients: 

_ / (mi +m 2 + h m s )n 

y min, m^n, . . . , m s n 

Fix a positive integer k > 2. We show that there exists a positive integer 
C(k) such that 



rin=i a fc« ^ 1 



■2 



for all positive integer t, if and only if GCD(m\, . . . , m a 



1. Mean Divisibility 



A sequence (a n ) (n = 1,2, ... ) of non zero integers is divisible if n \ m 
implies a n \ a m . It is strongly divisible if GCD(a n ,a m ) = \aGCD(n,m)\- Such 
divisibility attracts number theorists for a long time and a lot of papers dealt 
with properties of such sequences [131 EZ1 HI 13 [H E] . Especially primitive 
divisors of elliptic divisibility sequences and sequences arose in arithmetic 
dynamics are recently studied in detail OEUHj. In this paper, we introduce 
a weaker terminology which seems not studied before. We say that (a n ) is 
almost mean k-divisible, if there is a positive integer C = C(k) such that 
(lln=i a kn)/{Yli=i a n) £ for any positive integer t. In particular, (a n ) 
is mean divisible if IJn=i a " I \^i=i akn ^ or an y positive integer k and t. 
Clearly if (a n ) is divisible, then it is mean divisible. By definition, if a 
sequence is almost mean fc-divisible for all k with the constant C{k) = 1, 
then it is mean divisible. We are interested in giving non trivial examples of 
(almost) mean divisible sequences. In fact, we show that sequences defined 
by multinomial coefficients give such examples. Let mi, . . . ,m s be positive 
integers. A multinomial sequence is defined by 



a 



n 



(mi+m 2 H +m s )n\ ((mi + m 2 H hm s )n)! 



min, m2n, . . . , m s n J (min)\(m2n)\ . . . {m s n)\ 



Theorem 1. If GCD(mi,m2, ■ ■ ■ ,m s ) = 1, then the multinomial sequence 
is almost mean k-divisible for all k. 
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The proof relies on an interesting integral inequality (Lemma [3]) and its 
approximation by Riemann sums. Here are some illustrations: 

Corollary 2. 

Ilra=l ( 4n ) f- _}_n/ Y\n=l itn) a \n? [In=l (4n,8n,16n) ^ 

rr* r 5 ^ n rr' r 3 ^ ^ ' rr* r 7n \ 

\.i.n=l\2n) LLn=l\.n) lln=l U,2n,W 

/or any positive integer t. 

Readers will see that Figures 11121 and [3] in £JU] essentially tell why these are 
true. The constant C(k) is computed by an algorithm based on the proof 
of Theorem [TJ However it is not so easy to identify the set of t's at which 
the denominator actually appears. For the first example, there are infinitely 
many t with denominator 11, but the denominator 5 in the second example 
appears only when t = 2. See ^6]for details. We can also show 

Theorem 3. If GCD(mi,rri2, ■ ■ ■ ,m s ) > 1, then the multinomial sequence 
is not almost mean k-divisible for all k. 

Thus for a given k, a multinomial sequence is almost mean fc-divisible if 
and only if GCD(mi,m2, ■ ■ ■ , m s ) = 1 holds. For e.g., 

Corollary 4. Denominators of 

n* ( 8n ) 

LLn=l \4n) 
Iln=l Gn) 

for t = 1, 2, . . . forms an infinite set. 

The proof of Theorem [3] uses the Riemann sum approximation again, but 
we have to study more precisely the integral inequality of Lemma [3] at the 
place where it attains the equality. Indeed, we show that the set of primes 
in some arithmetic progression modulo k LCM{m\ , . . . , m s , Ylt=i m i) mus t 
appear in the denominators. If (a n ) is a multinomial sequence with param- 
eters (mi, . . . ,m s ), then so is (a£ n ) with (£mi, . . . ,£m s ). Since divisibility 
of (a n ) is hereditary to (ae n ), we see from Theorem [3l 

Corollary 5. Multinomial sequences are not divisible. 

Therefore multinomial sequences supply non trivial examples of almost 
mean /c-divisible sequences. For the mean divisibility, we can prove 

Theorem 6. If mi, . . . , m s are pairwise coprime and each nii divides Ya=i m i> 
then the multinomial sequence is mean divisible. 

For examples, we have 

Corollary 7. 

nt (2kn\ TT* ( 3fcn \ T-f* ( Gkn \ 

n=l \ kn) c ^ lln=l \kn,kn,kn) ^ lln=l \kn,2kn,3kn) ^ 

t 72nA~ G FT* ( 3n ~\1 J '( ^* 

n=l \ n) lln=l \n,n,n) lln=l \n,2n,3n) 

for any positive integer k and t. 
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The central binomial coefficient ( 2 ™) is of historical importance. By using 
( 2 n)' Chebyshev showed Bertran's postulate that there is a prime in any 
interval (n,2n]. Interesting divisibility problems on ( 2 ^) are discussed in 
|1CH [6j [11]. However the first example of Corollary [7J seems to be new. 
Theorem [6] models Theorem 5.2 in [2] which proves 

i-r2m (2j\ T-r2m+l (2j\ 

o-m 1 i j=m+l \ j / ^ o-m-1 Hj=m+1 V j / ™ 

Indeed the first example follows from Lemma 5.2 in [2] as well. 
Several further questions are exhibited in $7J 

2. Some Lemma 
Let mi, . . . , m s and fc be positive integers. Put m = and 



f(x) = [mx\ - ^[r 



\jriiX 

i=l 

We clearly have 

(1) f(x + l) = f(x). 

From x — 1 < |_^J < x, we see 



(2) [x\ + [~x\ 



-1 x 

iez 



Using ([2]), we obtain 

(3) /(x) + /(l-x) = a-l 

unless mx or mjX is an integer. Thus we have 



f(x)dx + / f{x)dx = {s- l)t. 
'l-t 



From this equality, 



/ (/(fee) - /(x)) dx = 0. 

./l-t 

Therefore we derive 

(4) / (f(kx) - f(x)) dx = 

J o 

and 

(5) f (f(kx) - f(x)) dx= I ' (f(kx) - f(x)) dx. 
Jo Jo 

First we assume s = 2. 
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Lemma 1. Let m\,rri2 be coprime integers and define 

f(x) = [(mi + m-2)xj — [w-ixj — |_ m 2xJ • 

Then for any positive integer k and any positive real t, we have 

r-t f-t 
f(kx)dx > / f(x)dx 
Jo 

and the equality holds if and only ift£ Ua^=o 



_ i a _i i 

k(mi+m,2) ' fc(mi+rrt2) 

Figure [T] and [2] in §[7] are the graphs of L* (f(kx) — f{x)) dx in special 
cases, which may help the reader. 

Proof. By ©, 

t 



(/(foe) - f(x))dx 







is invariant under £ i-> i + 1. We show the inequality for t € [0, 1]. The 
function f(x) is a right continuous step function with discontinuities at — Z 
and Ui=i with m = mi +m2. The discontinuities at ^Z gives +1 jump 
and (Ji=i gives —1 jump. Thus it is clear that /(x) = for x G [0, 1/m) 
and so / (x) = 1 for [(m — l)/m, 1) by ([3]) and right continuity. Since both 
integrands are identically zero, we have 

/ f(kx)dx = / f{x)dx 
Jo Jo 

for t £ [0, and the same is true for i G t^ 1 , 1] by ©. By (HJ, we first 
describe the distribution of the discontinuities in (0, 1), i.e., 

2 r ■ 
I i 

j = l,...,mi 



Recalling the idea of Farey fractions, let a/b and c/d be two non negative 
rational numbers with a, b,c,d£ Z, (a, b) = (c, <f) = 1 and a/b < c/d. Then 
(6) a a + c c 



b b + d d 



Arrange elements of 



2 , . 



j = l,...,mj -1 



in the increasing order. If j\/m\ and 32 /mi are adjacent, then we find 
a fraction {j\ + ji^/m in between. Considering the cardinality of -B, we 
notice that there exists exactly one element of B in the interval (i/m,(i + 
l)/m) for i = 1, 2, . . . , m — 2. For each i/m for i = 1, 2, . . . , m — 2, there 
is a right adjacent discontinuity of the form either j\/m\ or 32/^2- For 
the convenience, we formally extend this idea over the discontinuities at Z. 
There are no element of Ui=i ' m ~'^ J i n (( m — 1) / m ; 1) an d we associate ji/m\ 
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with j\ = mi for (m — l)/m and j%lmi with ji = m% for m/m = 1. In other 
words, we are formally treating in a way that 

m — 1 m\ 



m — 1 

m 



m m\ 



U 



m 777,2 
777 ' 7772 



though the last interval [m/m, 777-2/7772) is empty. We extend this convention 
to all positive integers by periodicity ([1]). Then for each discontinuity of the 
form — Z, there exists exactly one right adjacent discontinuity in U*=i ~'^ J - 
Thus we see that f(x) is a step function which takes exactly two values {0, 1}, 
and for any positive t, the integral f Q f(x)dx is computed as a sum of the 
length of intervals where f(x) = 1. These half open intervals have left end 

points in ^-Z and right end points in {t\ U [N?— 1 ^~ 
The required inequality is equivalent to 

ft rkt 

k / f(x)dx < / f(x)dx 



JO Jo 
and it suffices to show that 

G(t) = j f(x)dx -k f(x)dx 
Jo Jo 

is positive for < t < by fll]) and (j^J). Moreover the positivity 

is clearly true for -j^ < t < e and 1 — e < i < where e is the 

next discontinuity of f(kx) adjacent to -Ar. As /(a;) 6 {0,1} and G"(i) = 
- /(*))> w e see if f(t) = 1 then < and if f(t) = then 

G'(t) > 0. The minimum of Cr(£) is attained either at the end points of 
[e, 1 — e] or the point where G'(t) changes its parity from non positive to 
non negative, i.e., where / has negative jump. Thus it is enough to show 
that G{t) > for t 6 B, because the minimum must be equal to 

min{G(e), min G(t)}. 

t£zB 

Without loss of generality we may put t = j\/m\ with j\ = 1, 2, . . . , m\ — 1. 
Summing up the length of intervals where f(x) = 1, we have 

G -)= E-+E-- E - -ME-+E--E 

V 7771 / I ' 7771 7772 ^ — ' 777 / I 7771 7772 

with < j2 < 7772 and < •£ < fc — 1. The right side is 

(jf Jfe(fc - 1) + 2j' 2 M + £ 2 + e)ml- 2j 1 (j 2 (k - 1) + ^mim 2 + ^jfe(jfe - l)m| 

2777i7772(777l + 7772) 

whose numerator is a quadratic form of m\ and 7772 with the discriminant 

-4jjk£(k -i-l) 



^Here we think that the right adjacent discontinuity of a = — € — Z is the disconti- 
nuity ^ G -J-Z. 

J mo mo 
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and the coefficient of m\ is positive. Thus if < ^ < £; — 1 then G(^) > 0. 
For the remaining cases, we have 



2mim 2 (mi + m 2 )G 



_ \ k ( k ~ l)(j2"T-i - jim 2 ) 2 £ = 

mi) ~ I k(k - l)((j 2 + l)mi - jim 2 ) 2 £ = k-l 



Since mi and m 2 are coprime and j\ = 1, . . . , mi — 1, the right side can not 
vanish in both cases. We have shown the lemma. □ 

We prepare an elementary inequality: 

Lemma 2. Let s and b\, b 2 , ■ ■ ■ , b s +i be positive integers with GCD(bi, b 2 , . . . , b s 
1 and b = bi +b 2 + ■ ■ ■ + b s+ \. Then we have the following inequality: 



s+l / \ s+1 

11 GCD(b - bi, bi) GCD bj\ <J[ 

1=1 \ Kj<«+1 / 1=1 



b-h 



Proof. The left side is equal to: 

s+l 



S-tX I 

n 



i V 



GCD(b - b i+1 , bi+i) GCD b j 

l<j<3+l / 



with the convention b s+2 = b\. Then we see that GCD(b — 6j+i,6j+i) and 
GCD bj are coprime divisors of Therefore we have 

i<i<a+i 

s+l / \ s+l 

H GCD(b - b u h) CCD M < J] b k 

i=l \ i<j<s+i / k=l 

On the other hand 

s+l s+l / \ s+l 

i=l i=l \ j^i J k=l 

which proves the inequality. □ 

We wish to show a generalization of Lemma Q] 

Lemma 3. Let mi, . . . , m s be positive integers and put 

m = mi + m 2 + • • • + m s , g = GCD(mi,m 2 , . . . , m s ) 

and fix) = \mx\ — Yll=i l m i x \ ■ Then for any positive integer k and any 
positive real t, we have 



f{kx)dx > / f(x)dx 
Jo 

and the equality holds if and only if t £ Uo=o 



a 1 a I 1 

g km ' g km 
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Proof. Lemma [T] shows the case s = 2 and g = 1. The case s = 2 and 
5 > 1 is easily shown by applying Lemma Q] for m' = m/g and m- = rrii/g. 
We assume that the statement is valid until s(> 2) and prove the case 
m = mi + m 2 + • • • + m s+ \ and GCD(mi,m 2 , ■ ■ ■ , m 8 +i) = 1. The case 
GCD(mi,m 2 , • • • , m s _)_i) > 1 follows similarly to the case s = 2. 
By Lemma [2] with 6» = mj , we may assume 



(7) 



GCD(mi,m 2 , ■ ■ ■ ,m s ) GCD(y^^ rrij, m s+ i) < 



without loss of generality by changing indices. By the induction assumption, 
AO*) = L(E- = x m i) x \ — J2t=i l m i x \ satisfies the inequality: 



(8) 



fi(kx)dx > / fi(x)dx 



and the equality holds if and only if 



a=0 



a 

, — + 



9i k Ei=i m » ' 51 k E*=l 



m,- 



with 5i = GCD(mi,rri2, ■ ■ ■ ,m s ). Again by the induction assumption, for 
f 2 {x) = [mx\ - [{YH=i m i) x \ ~ [ m s+ix\ we nave 



(9) 



f 2 (kx)dx > / f 2 {x)dx 



and the equality holds if and only if 

b 1 



oo 

xe y 

6=0 



6 _ 1 

52 A;m ' g 2 km 



with 52 = GCD(^2 s i=1 mi,m s+ i). Since /(x) = f±(x) + /a (a;), we obtain 



f{kx)dx > I f(x)dx 



from (jHI and ([9]). Noting that 51,(72 are coprime, if either a/g± or 6/52 is 
not an integer, then \a/gi — b/g 2 \ > \/(gig 2 ). The inequality (Jj[J) shows that 



+ 



n 



b 16 1 

52 few ' 92 



7^ 



.51 ^Ei=i m *'5i k Yli=i m i. 
if and only if a/51 =b/g 2 G Z. We have shown the Lemma. □ 

Lemma 4. 27ie function f{x) in Lemma\M takes values in {0, 1, . . . , s — 1}. 

Proof. The case s = 2 is shown in the proof of Lemma [TJ Using the de- 
composition f(x) = fi(x) + /2(x) in the proof of LemmaU the assertion is 
easily shown by induction on s. □ 
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3. Proof of Theorem □ 

Let p be a prime and u p (n) be the largest exponent e such that p e divides 
n. Using the Legendre formula: u p (n\) = YlT=iY n /P e \-> we see * na t 

*(S^)-§£('®-'®) 

where / is defined in Lemma [3l Define 

t 



It suffices to prove that H(t) > for all t provided p e is sufficiently large. 
First we assume that p and k are coprime. Observe that 

t 



is a Riemann sum of the integral 



(f(kx)-f(x))dx. 

Our strategy is to show that these approximation is enough fine and Lemma 
[3] gives the answer to our problem. Since k and p are coprime, we see that 

(10) f(x) + /(l - x) = f(kx) + /(fc(l - x)) 

holds for x G prZ \ Z. Indeed by ([2D, both sides are equal to s — 1 + A(p e ) 
with 



A(p e ) = 5{p e | m) ~^5{p e 



1=1 

Here £(■) takes values 1 or according to whether the inside statement is 
true or not. We claim that 

(11) H{t) = H{p e -l-t). 

In fact, from (|10|) we see 

£('®-'®)^SL('®-'G 

and thus 

P e -i 



Therefore we have 

n=l x v ' v ' ' n=l 



p e J \ p e 



MEAN DIVISIBILITY OF MULTINOMIAL COEFFICIENTS 9 

which shows the claim. Put (x) = x — \x\. If (t/p e ) < 1/m, then H(t) > 
is clearly true because f(x) = for < x < 1/m. By (fTTl) . H(i) > also 
holds when {t/p e ) > 1 — 1/m. From Lemma there is a positive constant 
D such that 

rt 

(f(kx) - f(x))dx > D 

for 1/m < t < (m — l)/m. We have 

"t/ P e 

(f(kx) - f(x))dx 



Since f(kx) — f{x) is a step function, the last summand is zero if the interval 
((n — l)/p e ,n/p e ] contains no discontinuities, and its modulus is bounded 
from above by 2(s — l)/p e in light of Lemma [H Letting E be the number 
of discontinuities of f(kx) — f(x) in (0, 1), we have 



t/P \/(^-/(^-^E(/(^)-/(^ 

Therefore if p e > 2 ( S ~^) E then 



2(8 - 1)E 
< - 



2(s - 1)E 
> D - -i '— > 0. 

p e 



Exceptional discussion is required when p divides k. Putting k' = k/p Up ^ k \ 
we have 

£('®-'®)-S('(£)-'G 

since /(x) = for i£Z. If is a power of a prime p, then fc' = 1 and the 
right side is identically zero. If not, we have to replace k by k' and apply 
the same discussion. Then corresponding D' and E' are computed and we 
see that if K p := 2( " S ~^) E < p e then 



Therefore H(t) < happens only when 

2(s - 1)-B 

Z) ' p|fc and fc'>l 

which proves Theorem [U See ^6]for the actual computation of C(k). 



p e < M := max \ — , max K p , , 



10 



SHIGEKI AKIYAMA 



4. Proof of Theorem [3] 

We follow the same notation as in the previous section. It suffices to prove 
that there are infinitely many prime p's such that 

t oo 



S('l)-'(?))<" 



Take the minimum gap k between two adjacent discontinuities of f(kx). 
We shall find infinitely many such prime p's which are greater than 1/k and 
coprime with kLCM(m,mi,m2 ■ ■ ■ ,m s ). Since p > 1/k > km, a rational 
number with denominator p can not be a discontinuity of f(kx). The left 
side becomes 



£('(")-'©)■ 



if p > kt. By assumption g = GCD(mi,ni2, ■ ■ ■ , m s ) > 1, we know 

(f(kx)-f(x))dx = 

by Lemma [3j Hence from (|12p with e = 1 and 



for u G 

* = Lp/s. 



JL I + A. 

fern ' g km 



we have 



(13) 




fikx) 



/[- 



/(x) da; 



For the moment, we tentatively think that no discontinuities of fikx) inter- 
sects, i.e., 



3 

km 



j = l,...,fcm|(J(LI 



u=l 



j = l,...,kmi 



are the set of 'distinct' km + Yli=i = 2km points, and compute the right 
side. Since we are choosing a large p, there is at most one discontinuity £ in 
the interval ((n — l)/p,n/p\. If there is no discontinuity in ((n — l)/p,n/p], 

f(kx) = 0. If such a discontinuity £ exists, then we have 



then / ( f 



(14) 



±1 x€(^,0 



where ±1 is +1 if £ is the discontinuity of [fcmxj and —1 if £ is the discon- 
tinuity of \kmix\ for some i. Then we see 

(15) p (f(—)-f(kx))dx = ±(p^-(n-l)) = ±(pO. 



(n-l)/p 



P 
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A similar formula holds for f(n/p) — f(x). Summing up, from (Tl3j) we 
have shown: 



(16) 



km/g 

E 



P j \ 
km / 



/(- 

\P 



s krrii/g . 
1=1 j=l 




In reality, the discontinuities of f(kx) intersect in many places. For e.g., at 
least j/k for j 6 Z is a common discontinuity of L^w-i^J and L^w-j^J for 
i ^ j. However the above formula is correct without any changes. This 
is seen by a similar convention as in the proof of Lemma [U For e.g., if £ 
belongs to two discontinuities of \kmix\ and [kmjx\ with i ^ j, then 



'(t) 



/(fee) 



in ((n — l)/p, instead of (HH), but we computed integrand of (fT5|) twice in 

(Hi). 

Put = mi/g, m' = mj+m^- • -+m.g = m/g, L = LCM(m', m'i, . . . , m^) 
and L' = LCM(m'i, . . . , m^). An importance of the formula (I16p is that the 
value is determined by p (mod kLg). As p is coprime to kLg, by Dirichlet's 
theorem on primes in arithmetic progression, it suffices to show that there 
exists a single p such that the value of (fT6|) is negative to prove our theorem. 
Noting 



(17) 



km/g 

E 



/cm 



s krrii/g 

EE 

i=i j=i 




it is equivalent to show that there is such a p that 

fcm/g 



(18) E 



i=i 



PJ 
/cm 



s krrii/g 

EE 

i=i i=i 



PJ 

km; 



PJ 
m 



s rrii j g 

EE 

i=l j=l 



PJ_ 

mi 



is positiv^l. Moreover using ([2D, the value of (fTH|) changes its parity by the 
involution p -R- /cLg — p for p which is coprime to kLg, and the same holds 
for (|16p . Therefore our task is to show that either (|16|) or (|18p is not zero 
for some p which is coprime to kLg. We show this by dividing into three 
cases. 

Case m' > L' . One can find a p that 



P 



1 (mod kL g), p ^ 1 (mod /cm) 



2 From this expression, we see that (|16[) is integer valued. 
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and is coprime with kLg. From (|17p . the right side of <\16h becomes 

km/g . . . s krrii/g . (m/g . .. s rrii/g . 



^ \ km / EE fc m I E ( m ) EE 

i=l * ' i=l j=i % \j=i ' ' i=l j=i 

km/g , . . . . m/g 

P3\ 3 \ sr^ ( /P3\ 3 



rrij 



, . /cm / /cm / ^— ' V \ m / m 

j=i ^ / / j=i 



km/g 

(19) = E 

(mod fc) 



\ km I km 



km/g / " \ 1 fcm//ff 

\ I.- in / I.' i n — * ' 



Here we have 



\ /cm / km 
3 =l j=l 

j^O (mod k) j^O (mod k) 

where 7Z(a, b) is the minimum non zero integer congruent to a mod b. Be- 
cause the function (x/km) is increasing for < x < km, the minimum is 
attained by 

km/g km/g 

—— Y 7Z(j, km) = Y 

km ' ' /cm 

(mod fc) j^O (mod k) 

which shows that (|19|) is non negative. However, by p ^ 1 (mod /cm), clearly 
lZ(pj,km) takes values outside {lZ(j,km) \ 1< j < km, J ^ (mod /c)} 
and thus (|19|) must be positive. 
Case m' < L' . We choose a p that 

p ^ 1 (mod kL'g), p = 1 (mod /cm). 

From (fT7|) . the right side of (fTU|) becomes 

fcm/g . s tai/j . . . / m/g . s m^/g . 

j=i «=i j=i \j=i i=i j=i 

= -Y y (-)--)-yy(-)-- 

\\kmi/ kmij ^ ^ V\ m » / m * 

fcmi/g 

' ' P3 



E E 

i=i j=i 

j^O (mod k) 



kmt I kmi 



The inner sums are non negative and at least one of them is positive by the 
same discussion as in the former case. 
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CasE mf = V . We rewrite (|18p into 
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Hp) ■-- 



km/g 

E 



j^O (mod k) 



PJ 

km 



s km.il 'g 

E E 

i=i j=i 

(mod fc) 



PJ 

km; 



It suffices to show that there are two integers p\ and p2 which are coprime 
with km = k LCM(m, mi, . . . , m s ) and J~(pi) ^ J-{p2). 

First we study the case that there is a prime q with q | m' and q 2 \ km. 
Then we can take p\ = km /q — 1 and pi = km jq + 1 which are coprime with 
km. Since 



J | J 
q km 



J 

km 



j ^ (mod 

1 j = (mod 



we obtain, 



E i-E E 



l<j<km/g 
j^O (mod k) 
j=0 (mod q) 



l<j<km,i/g 
j^O (mod k) 
^-j=0 (mod,) 



Since m = m\ + 777-2 + • • • + m s and m/mi are integers, the right side is 
expected to be negative. However a careful computation is required, because 
not all kmi/g are divisible by q. Since q divides either m\ or m/mj = m' /m' i7 
we have 



E 



l<j<krrii/g 
j^O (mod k) 
^j=0 (mod?) 



km\ — m\ 

km'i km\ 
y ~g LCM(k,q) 



m/m-i =0 (mod q) 
m/mi ^ (mod q) 



3 Thereare such pairs, for e.g., 1 + 2 + 3 = LCM(1,2,3), 2 + 3 + 3 + 4 = LCM(2, 3, 3, 4). 
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Therefore 

F(P2) -F{pi) 



km 1 


km' 


Q 


LCM{k,q) 


km'i 


km\ 


. 1 


LCM(k,q)_ 


km\ 


km\ ^ 


Q 


LCM(k,q) y 



Ki<s >^// l<j<s 

^0" (mod g) -2- =0~ (mod g) 



l<i<s 
2-=0 (mod g) 



(A; - l)((LCM(fc, c?) + - 1) + 1) - LCM(k, q) + q 



i£ ^ *LCM( M ) 

^=0 (mod,) 

< V (g + 2 ) m » < n 

^ qLCM(k,q) ■ 

^=0 (mod?) 

Here we used k>2 and q > 2. Note that the last sum is non empty, because 
q | m' and GCD(m' 1 , . . . , m' s ) = 1 implies that g divides m/m, = m'/m^ for 
at least one i. 

Second, consider the case that there is a prime divisor g > 3 of m'. We 
may assume that either 

/ km km \ ( 2km 2km 

(Pi,P2) = 1, + 1 or 1, + 1 

V q q ) V Q Q 

satisfies GCD(pi, km) = 1 for i = 1 and 2. In fact, if for e.g. km/q —1 = 
(mod q) and 2km jq + 1 = (mod q) hold, then q 2 divides km, which is 
reduced to the first case. Once we have such a pair (pi,P2), we can show 

Hpi) - Hpi) < o 

in the same manner. So we finally consider the case that all the prime 
divisors of m! is 2 and 3 and m' is square free, which covers the remaining 
cases. There exists only one such case with 

m'j H h m' s = LCMim'i, m' s ), 

that is, s = 3 and (m^m^m^) = (1,2,3). So our last task is to consider 
the case: (mi, m.2, m$) = (g, 2g, 3g). Since m = 6g is even we can choose: 



(Pl,P2) 



- 1, *p + 1) fcm = (mod 4) 



(*P - 2, *p + 2) fcm = 2 (mod 4) 
Then we see F(p2) — F{pi) < in the same manner. 
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5. Proof of Theorem [6] 

This proof is inspired by Theorem 5.2 and Lemma 5.2 in [2]. We use 
the same terminology as in section [3l Under the assumption, f(x) has no 
discontinuity of negative jump in (0, 1) by cancellation. Thus f(x) is non 
decreasing. We show that 

where k' = k/p Up ^ k \ By periodicity of /, it suffices to show the case that 
t < p e . Since k' and p are coprime, k'n mod p e (n = 1, 2, . . . ,p e — 1) are 
distinct. From f(l/p e ) < f{2/p e ) < < f{{p e - l)/p e ), we see 




is the minimum of the sum of t elements in {f(i/p e ) \ i = 1, 2, . . . ,p e — 1}, 
which finishes the proof. 



6. Computation of C(k) 

In this section, we explain the computation of the smallest C(k) by an 
algorithm based on the proof of Theorem [TJ For a given k, first we compute 
the minimum D q of I(t) = f{k'x) — f{x)dx for t G [1/m, 1 — 1/m) for 
k' = k'(q) = k/q u vW for all prime divisor qoik. Denote by D\ the minimum 
of I {f) for k' = k. Since the minimum of I(t) is attained at the discontinuity 
of the step function: I'{t) = f(k't)—f(t), using (fTTj) . it is explicitly computed 
as 



minl(t) = min < 



min I — 



mm mm 1 I — — 

Ki<s k f m^ ^ k'rrij \ fx, Tflj 



By Theorem[Tl we know D q > 0. Number of discontinuities is bounded from 
above by 2&/771Q. Then we compute for prime p's with 

(2km(s-l) 2k'(q)m(s-l) 
p < M := max ^ ) -, max — — 

[ U\ q\k and k'>l Dq 



the values 



r p h_i * h+e-i 

fi(p) = min < 0, min 

I ~ n=l e=l 



k'n 

p e 




It is better to take the exact value to make faster the computation, as we do below in 



examples. 
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with h = [_log(M)/log(p)J and £ = |_log(fc'm)/log(p)J . By the proof of 
Theorem [H if p e > M then 

±0(2) -t 



n=l 



> 



k'n ^ J_ 



for any positive integer t. Further by ^^pj 



£('G?)-'te 

for e > /i + £ and i < p . We have 

t oo 



fx(p) = min jo, mm £j (/ (^r) " / {je 

and attained as the denominator for some t. We obtained the best 

constant C(Jfe) = n p <Jif P** 60 . D 

We briefly demonstrate this algorithm by showing Corollary [2] and make 
precise the comments afterwords. For s = 2, mi = 2, m,2 = 3 and k = 2, 
the graph of the function for k' = k \s depicted in Figure [TJ We may 



0.02 



Figure 1. mi = 2, m 2 = 3, k = 2 

take .Di = 1/60 and £i = 8. Checking all primes which do not exceed 
2 • 8 • 60 = 960, we found the only non zero output = —1. We can 

confirm that 

t oo 3 

£ £ (/(2n/ll e ) - /(n/ll e )) = £ (/(2n/ll) - /(n/11)) = -1 

n=l e=l n=l 

when i = 3 + ll fc and k > 2, since /(x) = for < x < 1/5. Therefore the 
denominator 11 in the first formula in Corollary [2] appears infinitely often. 
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For s = 2, mi = 1, m,2 = 2 and k = 3, the graph of I(t) is depicted in 
Figure [2j We have D\ = 1/18, Si = 10 and the only non zero //-value is 



0.15 - 



0.10 - 



0.05 - 




0.2 0.4 0.6 0.8 1.0 

Figure 2. m\ = 1, m 2 = 2, = 3 
/i(5) = — 1. In this case, one can also confirm that 

2+11* f 

^(/(3n/5 e )-/(n/5 e )) = -l. 

n=l e=l 

However we have 

i /+i 

^^(/(3n/5 e )-/(n/5 e ))>0 

n=l e=l 

for 2 < t < 5* — 1 and any positive integer /. In other words, the function 
Sn=i Se=i (/(3 ra /5 e ) — /(n/5 e )) has a period 5^ and attains —1 infinitely 
often as above, but such negative values are erased by the next period of 
length 5^ +1 , except when t = 2. The denominator 5 in the second formula 
of Corollary [2] appears only when t = 2. 

For s = 3, mi = 1, m2 = 2, 7713 = 4 and fc = 4, the graph of is depicted 
in Figure O We have D\ = 11/112, i?i =32 and fj,(p) = for all prime p. 

In this manner, a prime divisor p of C{k) actually appears in the denom- 
inator only when the p-adic expansion of t has a special form, and not easy 
to describe the set of such t's. 

7. Questions 

We wish to list several open problems. A sequence (a n ) may be called 
almost mean divisible if there is a positive integer C such that 

t t 

(II a ^)/(II a ™) e c Z 

n=l i=l 
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Figure 3. mi = 1, mi = 2, m% = 4, k = 4 

for any positive integers k and i. In other words, (a n ) is almost mean divis- 
ible, if it is almost mean fc-divisible with a uniform constant C independent 
of the choice of k. 

• Is there an almost mean divisible multinomial sequence, which is not 
mean divisible ? 

• Is there a mean divisible multinomial sequence which does not satisfy 
the condition of Theorem [6] ? 

• Is there any other (almost) mean divisible sequence of number the- 
oretical interests ? 

I expect the answer for the first question is negative, because the bound 
M increases as k becomes large, in the proof of Theorem [TJ For the second, 
there may exist such multinomial sequences for s > 2. For s = 3 and 
(mi, 7712,1713) = (m, 1, 1), I checked by the algorithm in $6] to obtain a 

Corollary 8. 7/3 < m < 10 and k < 10, then 

t t (m+2)fcn \ 

n\kmn,kn,kn) nj 
((m+2)n\ G 
n=l \mn,n,nJ 

hold^ for all positive integer t. 

We do not know if this is true for all k for some m > 3. This sequence is 
factored into two: 

\m + 2)n\ _ ((m + 2)n\ /2ra N 
mn,n,n J \ 2n J \n 

The former sequence (( (m + 2)n )) is almost mean fc-divisible for all k by The- 
orem [Tj and (( 2 ™)) is mean divisible by Theorem [6j So the denominators 
generated by the first sequence might be canceled by the numerators from 



^By Theorem [6] this is valid for all positive integer k when m = 1 and 2. 
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the later one. It is an interesting problem to characterize all mean divisible 
multinomial sequences. 

As for the third question, we can construct a different type of non divisible 
almost mean fc-divisible sequences. Fix an integer I > 1 and let a and /3 be 
conjugate quadratic integers so that a/j3 is not a root of unity. Define the 
£-th homogeneous cyclotomic polynomial: 



*t(x,y)= [] {x-Cv) 



0<m<£ 
GCD(m,£)=l 

where £ is the primitive £-th root of unity. Put 

Then (c n ) is a non zero integer sequence and for any integer k coprime to 
£, we have c n | c kn . This implies that (c„) is almost mean /c-divisible for 
GCD(kJ) = 1. For example, taking I = 2, 

^(¥)" + (¥)" 

gives the Lucas sequence, which is non divisible but almost mean /c-divisible 
for all odd integer k. However, it may not be a significant construction 
because they already have divisibility c n \ c kn not for all but for some k. 

Acknowledgment. The author would like to express his gratitude to 
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